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I. INTRODUCTION 

Establishing the thermodynamic cost of various op- 
erations processing information is a fundamental issue 
technologically as well as conceptually. Its origin can 
probably be traced back to Maxwell discussing his well 
known demon [1 . An important result in this respect was 
achieved by Landauer who established a lower bound for 
the heat generated when erasing one bit of information 
[2]. In doing so, he established an explicit connection be- 
tween a thermodynamic quantity, heat, and a quantity 
measuring information, the Shannon entropy [3 j. The va- 
lidity of this so-called Landauer principle has been ver- 
ified in the context of classical, quantum and stochas- 
tic dynamics [UJS]. Its relation to entropy production 
and microscopic reversibility was discussed in Refs.|7H9]. 
The implications of Landauer's principle for computation 
were recognized early on [T0HT2] and opened the way to 
the field of reversible computing [13] . 

With the advent of stochastic thermodynamics [HJ[T5], 
Landauer's principle has become an immediate conse- 
quence of the second law. In this formalism, a sys- 
tem consists of states i with energies and probabil- 
ities pi. The Shannon entropy of the system is given 
by S = — kB^ZiPi^Pi- When the system is in con- 
tact with a single reservoir at temperature T, transi- 
tions between the system states occur and the system 
probabilities evolve according to the Markovian master 
equation pi = . WijPj with transition rates Wij which 
satisfy the local detailed balance condition Wij/wji = 
exp {(ej — €{)/ (k B T)}. The second law in stochastic ther- 
modynamics reads AS = Q/T + A[S, where AS is the 
change in Shannon entropy, Q is the integrated heat flow 
Q = €ifti entering the system, and A\S is the nonneg- 
ative entropy production which only vanishes for quasi- 
static transformations where detailed balance is satisfied. 
If the system is a bit (i.e. a two level system with two 
states and 1) initially containing the maximal infor- 
mation S = k B \n2 corresponding to the uniform prob- 
ability po = pi = 1/2, Landauer's principle states that 
erasing that information, i.e. bringing the initial system 
entropy to S = 0, will produce an amount of heat of at 
least k B T hi2. This immediately follows from the sec- 
ond law since AS = k B In 2 and therefore the generated 
heat reads -Q = k B T ln2 + A { S > k B T ln2. Landauer's 
lower bound is only reached for quasi-static transforma- 



tions where A\S = and thus requires an infinite amount 
of time. 

Since stochastic thermodynamics naturally combines 
dynamics with thermodynamics, it opens the way to the 
study of information erasure in finite-time. Interesting 
results in this direction have been obtained for systems 
described by Fokker-Plank equations. Transformations 
of duration t between two sets of probabilities which min- 
imize the heat generated lead to an entropy production 
scaling as 1/t [T6l [17]. This result also holds for sys- 
tems described by master equations but is limited to a 
regime of low dissipation [18] . These studies have impor- 
tant implications for the study of efficiencies in finite-time 
thermodynamics [16, 19-22 . Furthermore, many recent 
works have analyzed the implications that feedback con- 
trol may have on the thermodynamic description of a 
system [23H32] . In this paper we are going to build on 
these studies to analyze the process of information era- 
sure in finite time first without and then in the presence 
of feedback control. 

In the first part of this paper, i.e. section [il] we intro- 
duce our model and study in detail the erasure of infor- 
mation in finite time by analyzing the trade-offs between 
generated heat, erasure time, and accuracy of the era- 
sure. We also introduce the notion of erasing efficiency 
and erasing power. In the second part of the paper, i.e. 
section |III[ we study how to improve the erasure pro- 
cess by introducing a feedback process. Conclusions are 
drawn in section HVl 

II. ERASING INFORMATION IN FINITE TIME 

We consider classical information stored in an array of 
single level quantum dots. Each dot constitutes a classi- 
cal bit since it can either be empty or filled with an elec- 
tron (0 or 1) with probability 1 — p and p respectively. 
The Shannon entropy per bit is 

S = — fcsplnp — k B (l — p) ln(l — p). (1) 

It takes its maximal value S = k B In 2 when p =1/2 and 
its minimal value S = when p = or p = 1. The energy 
of an electron in the dot is denoted by E. The stored 
information is metastable in the sense that the energy 
gap E — fi env to bring an electron in or out of the dot from 
or in its surrounding environment is much larger than the 
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energy fluctuations k B T of the environment. However, it 
can be modified when the dot enters in contact with a 
metallic lead moving at constant speed along an array of 
quantum dots as depicted in Fig. [I] The speed of the lead 
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FIG. 1. (Color online) a) The information stored in the array 
of single level quantum dots is erased by putting the quantum 
dots one after another in contact with the metallic lead and 
applying, during the contact time t, a time dependent proto- 
col on the lead chemical potential /i(t). b) The decrease of 
the initial Shannon entropy Si = ks log 2 is accompanied by 
an heat release — Q in the environment. 

controls the contact time t between the dot and the lead. 
The lead is at the surrounding temperature T but its 
chemical potential ji(t) is externally controlled. Its time 
dependence during the contact time t will be denoted as 
the protocol and is the same for each dots. The dynamics 
of the dot during the contact time is described by the 
master equation 



p=-C(l-f(E))p + Cf(E)(l-p), 



(2) 



where Cf(E) (resp. (7(1 — f(E))) is the rate at which 
the lead can donate (resp. receive) an electron to (resp. 
from) the dot, and f(E) = (exp {{E - n)/(k B T)} + l)" 1 
is the Fermi distribution of the lead. Introducing the 
variable e = E — ji, setting fc# = 1, and measuring time 
in units of C _1 , we can rewrite Q as 

^-^exp^vn+i' (3) 

When solving this equation over the contact time t with 
a time dependent protocol e(t), the probability p evolves 
from an initial value Pi = p(0) to a final value Pf = p(t). 
The resulting change in Shannon entropy per bit is given 
by the second law of stochastic thermodynamics [33] 



AS = Sf - Si = ^ + AiS, 



(4) 



where Sf and Si are the Shannon entropies corresponding 
to pf and Pi respectively. The heat entering the dot is 
given by 



dr p(r)e(r) 



(5) 







and the resulting nonzero entropy production reads 
'•' 1-P(r) e(r) 



A iS = J o rf T p( T )(lni-^-^)>0. (6) 



An erasure process is characterized by a negative en- 
tropy change AS < 0. Due to Eq. Q, this process re- 
leases heat into the environment (i.e. the heat absorbed 
is negative Q < 0). We define the efficiency of this pro- 
cess as the amount of entropy change resulting from this 
heat release 



< 77 
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< 1. 



(7) 



The upper bound of the efficiency 77 — 1 corresponds to 
Landauer's lower bound and is reached when the proto- 
col evolves quasistatically between e(0) = T ln(l /pi — 1) 
and e(t) = T\n(l/pf — 1) so that the entropy production 
vanish Ai*S = 0. In this case the erasure process becomes 
infinitely slow. The lower bound rj = will be reached 
when the heat generated for a given erasure diverges. 
Since the second term in Eq. ([3| is bounded between 

and 1, the fastest way to decrease (resp. increase) the 
probability p is given by p(t) = —p(t) (resp. p(t) = 

1 —p(t))- This corresponds to a protocol e — > oo (resp. 
e — » — oo ) which leads to a divergent heat production and 
thus a vanishing efficiency 77 = 0. This argument also 
implies that the time required to bring the probability 
from pi = p(0) to pf = p(t) is always larger or equal to a 
minimal time tm\n- 



if pi > pf : t > t n 
if pf > pi'. t>t n 



In 



In 



Pf' 

1 — Pi 
1-Pf' 



(8) 



The case t = t m i n corresponds to the protocol e — >• 00 
which leads to 77 = 0. If t < t m i n , no protocol is able to 
reach pf from p it This brings us to the important result 
that perfect erasure (Sf = due to p/ = or = 1) 
leads to a divergent t m i n . In other words it is impossible 
to completely erase an initial finite Shannon entropy per 
bit in a finite amount of time. Only quasistatic process 
can do so, since in this case the contact time and the 
minimal time can diverge preserving t > t m i n . 

We can of course reformulate (|8| by saying that for a 
finite contact time t, the final probability pf which de- 
fines the erasure error Sf is always bounded by a critical 
probability p c : 

if Pi>p f : Pf >p c = Pi exp {-t}, 

if pf > pi : pf < p c = 1 - (1 - p^ exp {-£}. (9) 

As a result, minimizing the erasing error for a given con- 
tact time corresponds to the case Pf = p c which leads to 
a divergent heat and thus to a vanishing efficiency 77 = 0. 

We now consider finite contact times t with fixed Pi and 
Pf such that Pf < p c and try to find the protocol which 
maximizes the efficiency. Since AS is fixed, the highest 
efficiency will be obtained when the heat released in the 
environment is minimal or equivalently when the entropy 
production is minimal. The procedure to find the optimal 
protocol minimizing the heat released is detailed in Ref. 
[T8] . The amount of heat generated with the optimal 
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FIG. 2. (Color online) Minimal heat generated —Q/T when 
erasing in a finite time t an initial information Si — Tin 2 
(jpi — 1/2) with a remaining error Sf (pf). The hori- 
zontal (red) curves on the surface correspond to Q/T = 
0.05, 0.3, log 2, 0.8. The vertical plane (orange) corresponds 
to pf = p c where the heat diverges. The (black) curve in the 
t = 20 plane corresponds to Landauer's lower bound where 
Q/T = AS. 



protocol bringing the initial probability pi = 1/2 to the 
final value pf (corresponding to an erasing error Sf) in a 
time t is displayed in Fig. [2] The corresponding erasing 
efficiency is depicted Fig. j3] As the contact time in- 
creases and the protocol approaches the quasistatic solu- 
tion, the heat approaches the Landauer limit Q/T = AS 
represented by the full (black) line in the t = 20 plane of 
Fig. [2] and the erasure efficiency increases and approaches 
one in Fig. [3j Also, as the final probability pf reaches 
its critical value p c , the heat starts diverging and the 
efficiency drops to zero. This region where pf — p c is 
represented in Fig. ^1 by a vertical (orange) plane and in 
Fig. [3] by the full (black) line in the 77 = 0.2 plane. In 
the limit pf — >• 1/2, not surprisingly the heat vanishes, 




FIG. 3. (Color online) Erasing efficiency 77 corresponding to 
the erasure process in Fig. [2] The horizontal (red) curves on 
the surface correspond to 77 = 0.3, 0.5, 0.7. 



but the efficiency converges to the full (black) line in the 
Pf = 0.5 plane of Fig. 13] which can be calculated analyti- 
cally using results of Ref. [18] and gives 77 = (1 + 2/t) — x . 
In Fig. |2j the behavior of pf as a function of t for a 
fixed value of heat Q is shown by horizontal (red) curves. 
The long time behavior of these curves depends on the 
value of the heat. For — Q < Tin 2, it decreases to an 
asymptotic value of pf corresponding to the Landauer 
limit AS = S f - In 2 = Q/T. However, for -Q > Tin 2, 
it will eventually reach pf = 0. 

In order to quantify the average amount of informa- 
tion erased per unit time during a contact time t with a 
given amount of generated heat — Q, we define the era- 
sure power 

— A 9 

r(Q,t) = — , (10) 

which is a function of Q and t. The Landauer limit which 
leads to an optimal erasure efficiency (77 = 1) corresponds 
to zero erasure power (V = 0). Nonzero erasure power 
only occur at finite contact time as shown in Fig.[4]where 
V is plotted as a function of the heat generated (with the 
optimal protocol minimizing heat) and time. We note 




FIG. 4. (Color online) Erasure power V as a function of heat 
generated —Q/T and the duration of the erasure t. The (red) 
curves on the surface correspond to Q/T = 0.2, 0.7, 1.2. 

that for a constant generated heat —Q/T, the erasure 
power V reaches a maximum value for relatively short 
times and then drops to zero in the long time limit where 
the optimal protocol becomes quasistatic. 

The erasure efficiency at maximum erasure power 77*, 
as well as the corresponding contact time £*, the corre- 
sponding final value of the probability p £ , and the corre- 
sponding value of the power V* are displayed in Fig. [5) 
In the limit of small generated heat (—Q/T — » 0) the en- 
tropy change has to vanish since the second law imposes 
— AS < —Q/T. It does so as (1/2 - Pf) 2 since in that 
limit the final probability pf approaches 1/2 which is a 
maximum of AS. In the same limit the contact time can 
be shown to behave as 1/2 — pf so that the erasure power 
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FIG. 5. (Color online) From upper to lower frame: Erasing 
efficiency at maximum erasing power (black) , the correspond- 
ing contact time (green), final probability (red), and power 
(blue). 



and the corresponding efficiency both vanish. In the op- 
posite limit of large generated heat (— Q/T — » oo), using 
(|8|, the erasure power behaves as V = — AS / hi(pi / p f) 
and its maximum value occurs at pf « 0.10892. As a con- 
sequence the efficiency decreases in that limit as r] oc 



III. ERASURE WITH FEEDBACK 

We now turn to an erasure process assisted by the feed- 
back process depicted in Fig. [6] An imperfect measure- 
ment is performed on the bit to be erased and the ensuing 
protocol depends on the output of that measurement. 
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FIG. 6. (Color online) a) As the lead enters in contact with 
the quantum dot to be erased, an imperfect measurement 
is performed, b) The time dependent protocol p(t) applied 
during the contact time t depends on the output of that mea- 
surement. 



The two possible states of the bit, empty and filled, are 
denoted by a = 0, 1 and the probability to find the bit in 
state <j at time t is denoted P t (a) . We consider ideal mea- 
surements which do not perturb the system measured. 
Therefore, the probability to find the bit in a given state 
a remains the same right after as right before the mea- 
surement which occurs at time t = 0: Po(a = 1) = pi 
and Po(<j = 0) = 1 — Pi. The two possible outcomes of 
the measurement are denoted by a = 0, 1. The accuracy 
of the measurement apparatus is characterized by the 
conditional probability Po(a\a) to find the bit in state a 
when the measurement outcome a is realized: 



1 — S if a = a 
S if a ^ a 



(ii) 



A perfect measurement corresponds to S = and fully 
characterizes the system state while a useless measure- 
ment corresponds to S = 1/2 and does not provide any 
additional information about the system state. After 
the measurement, the protocol e(<r) depends on the mea- 
surement outcomes. This means that the system will be 
described by the conditional probabilities P t (a\a) which 
evolve according to the master equation Q with the ini- 
tial condition (11). The joint probability P t (a,a) is re- 
lated to the time-dependent conditional probability by 



P t (a,a) = P t (a\a)P(a). 
It corresponds to the Shannon entropy 

S t (a, a) = -J2 P t a) ^ Pt {?, a) 



(12) 



where in the second line we defined 

S(<x) = -£P(<x)lnP(<7) 



S t (a\a) = -^2P t (a\a)lnP t (a\a). 



(13) 

(14) 
(15) 



The probability to measure an output cr, P(cr) in (12), is 
obtained from the condition 



Po(a) = Y^P Q {a,a) = £p (<t|<t)P(<t). 

a a 

We easily find that for our model 

P(a=l) = l-P( ( ? = 0) = ^. 
The Shannon entropy of the bit at time t, 
S t (<x) = -^P t (<x)lnP t (<x), 



(16) 



(17) 



(18) 



is related to the mutual information between the system 
and the measurement outcome by [34] 



M t = S(a) + S t (a)-S t (a,a) 
Pt(v,<r) 



(19) 



Pt(v)P(v) 



> 0. 



5 



By combining Eqs. (13) and (19), the change in the bit 



entropy can be written as 

AS t = AM t + ^P((t)ASt((T\<t). (20) 

Since P t (a\a) evolves according to Eq. ^ with protocol 
e(a), using traditional stochastic thermodynamics at the 
level of this conditional probability, we find that 



AS t (a\a)=Q(a)/T + A i S(a), 



(21) 



where Q(a) and A[S(a) > are the heat and the en- 
tropy production associated to the dynamics following a 
measurement output a. We can thus rewrite (20) as 



AS t = AM t + Qf/T + (AiS), (22) 
where AM t is the change in mutual information and 

Q F = Y^P{cr)Q{a) (23) 

{AiS) = ^2 P (a) AiS (a) > 

a 

are respectively the heat and entropy production aver- 
aged over the possible measurement output giving rise to 
different protocols. Eq. (22) can be seen as a general- 



ization of the second law of stochastic thermodynamics 
Q in absence of feedback to situations with feedback. It 
imposes the following bound on the heat released by the 
system 



Qf > TAM t - TAS t > -QF,min, 



(24) 



where the minimum heat release possible is given by 



Qf,. 



-TMq - TAS t . 



(25) 



In order for the released heat —Qf to reach the interme- 



diate bound in (24), all the erasure processes following 



the measurement have to be performed quasistatically: 
AiS {a) = for all a. The bound —Qf, m i n can only 
be reached if in addition all the different protocols end 
up at time t at a same final value independently on the 
measurement outputs. Indeed, since the probability of 
a quasistatic processes is fully determined by its proto- 
col, at the end of the each process we would have that 
P t (a\a) = Pt(cr) : and thus that all the initial mutual in- 
formation has been consumed at the end of the process: 
M t = 0. 

We will assume, as we did in last section, that pi = 1/2. 
Therefore, using (17), we have P(a = 1) = 1 — P(a = 0) 
= 1/2. With the measurement errors (11) and using (17) 
and (fl9|, the initial mutual information becomes 



M = In 2 - S 6 , 

Ss = SlnS-(l-S) ln(l-S). 



(26) 



We have seen in last section that perfect erasure can 
be achieved in the quasistatic limit. This means that 
for each of the two measurement output, we can reach 



Pt(a = l\a) P t (a = 1) = so that AS t = - In 2 and 
M t = 0. As a result, the minimum possible heat released 
in presence of feedback is given by 



?F,min = Tin 2 — TMq = TS$. 



(27) 



This result can be viewed as an extension of the Landauer 
principle in presence of a feedback. In the limit of a 
perfect measurement (5 = 0) the minimal heat released 
completely vanishes. 

We turn now to the protocols €t(cr) that minimize 
the released heat required to erase a given amount of 
information in finite time. To make the comparison 
with last section meaningful, we impose the same change 
in Shannon entropy with and without feedback. We 
will therefore minimize the heat released in going from 
Po(cr = 1) = Pi = 1/2 to = 1) = pf in a finite 

time t. The final values of the conditional probabilities 
Pt(<J — have to satisfy the constraint 



P t (a) = Y J PM°)P{°)- 



(28) 



To simplify the notation we define 

pf =P t (a = l\a) = l-P t (a = 0\a), (29) 
Since P{a) = 1/2, we find that 

P f ) =2p f -pf. (30) 

This fixes p^ in terms of pf but leaves p^p free. The 
average heat released now reads 

Q F = Y^P{a)Q{a)= 1 -{Q{Q) + Q{\)). (31) 

(J 

Minimizing —Qf for a given time t is done by first sep- 
arately finding the two optimal protocols e t (0) and e^(l) 
minimizing Q(0) and Q(l) respectively, using the final 
probabilities p^ and p^ which are related by (30). The 
second step consists in further minimizing the resulting 
expelled heat —Qf with respect to p^\ The minimum 

is reached for pf\ pt • 

In absence of feedback we have seen that to avoid di- 
vergences in the heat we need to fulfill the condition 
2pj > exp{— t}. Similarly in presence of feedback, to 
avoid divergences in the heat Q(0) and Q(l), we need 
to satisfy p^ > (1 — S)ex.p{—t} and p^ > Sex.p{— t} 
respectively. Using p0| ), these two conditions combine as 

(l-£)exp{-£} <p { p < 2p f -Sexp{-t}. (32) 

In the limit 2pf — » exp {— t} where the heat of the process 
without feedback becomes divergent, the bounds collapse 
and all the critical final probabilities become identical. 
For 2pf > exp{— £}, P^ opt will be located within the 
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ing (21) with (15), we get that 




FIG. 7. (Color online) Ratio Qf/Q as a function of the fi- 
nal probability pf for different erasure times t. The curves 
bend down as t increases and the dashed curve corresponds 
to the quasistatic limit (35). The vertical lines denote the 
critical probability p c reachable for each erasure time t 
measurement error is fixed at 5 — 0.02. 



The 



bounds ( 32 ) where none of the critical final probabilities 



are reached. 

In the quasistatic limit, we can derive the minimal 
value of —Qf predicted by (25) with (26) explicitly. Us- 




FIG. 8. (Color online) Ratio Qf/Q as a function of the mea- 
surement error 5 for different erasure times t. The final prob- 
ability is fixed at pf = 0.01. The dotted curve represents the 
ratio at the critical time t m i n = —\n2pf while the dashed 
curve corresponds to the quasistatic limit (35). 



Qf = \ {-pf Inpf - (1 - pf) ln(l - pf) - S S ) 

+ l(-pfln P f-(l- P f)Ml-pf)-S s ). 

(33) 



If we impose the constraint (|30|), the derivative of 

jo) 

*7,opt 

AMt-AS u 



with respect to vanishes for P^ opt 



Pf- 

Since in the quasistatic limit, —Qf/T = AM t — ASt, by 
minimizing the heat at fixed AS, we are in fact minimiz- 
ing the final mutual information M t which vanish pre- 
cisely when p^p = p^p . 
of — Qf then reads 



Using (26), the minimum value 



>F,min = —TM - TAS t = TS S - TS 



(34) 



A relevant quantity to look at is the ratio Qf/Q be- 
tween the heat released in presence and in absence of the 
feedback to erase a given amount of information in finite 
time. In the quasistatic limit, using (34), we find that 



Qf,v 



Qn 



m2-sy 



(35) 



Since S$ < In 2, except for useless measurements (S = 
1/2), the feedback always helps to reduce the expelled 
heat. In finite time, the ratio Qf/Q corresponding to the 
protocols minimizing heat has been calculated in Fig. [7| 
for a given measurement accuracy, as a function of the 
final probability pf and of the erasing time t. This ratio 
is always lower than one and decreases for longer erasure 
times. This shows that the reduction in the expelled heat 
thanks to the feedback increases as the contact time in- 
creases and is most significant in the quasistatic limit. 
The same ratio has been calculated in Fig. [8j for a fixed 
value of the final probability pf = 0.01, as a function of 
the measurement error S and of the erasing time t. As 
we approach the useless measurements limit 5 — >• 1/2, 
the ratio goes to one independently of the erasure time 
indicating no gain by the feedback. In the perfect mea- 
surement limit S — >• 0, as the contact time becomes longer 
the ratio approaches zero and may even become negative. 
This happens in the t — > oo limit when Ss < Sf as pre- 
dicted by 434b. 



IV. CONCLUSIONS 

We proposed in this paper a model to study the ther- 
modynamics of information erasure in finite time. It con- 
sists of an array of single level quantum dots which can 
store classical information given that each dot with its 
two states, with or without an electron, constitutes a 
classical bit. The initial information is quantified by the 
Shannon entropy Si. The erasure process consists in de- 
creasing that information by — AS = Si — Sf > and is 
performed by a metallic lead moving at constant speed 
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along the array of quantum dots. During each lead-dot 
contact time t a time dependent protocol controls the 
lead chemical potential. 

In the first part of the paper we considered the situa- 
tion where the same protocol is applied to each dot. We 
found the following results. Perfect erasure Sf = of 
any finite initial information Si always requires an in- 
finitely long contact time t — >• oo. Imperfect erasure 
(which leaves some errors Sf) with a maximal erasing 
efficiency 77 = 1 requires a quasistatic protocol and thus 
again a diverging contact time. The erasing efficiency 
r] measures the fraction of the generated heat (divided 
by temperature) (-Q/T) that is used to erase the in- 
formation while the rest is lost as entropy production. 
Reaching Landauer's lower bound Si = —Q/T requires 
perfect erasure with efficiency one. In finite time, errors 
and entropy production are unavoidable. Attempting to 
minimize errors will result in a diverging heat and a van- 
ishing erasing efficiency. For larger values of the error, we 
calculated the protocols minimizing the heat generation 
and studied the resulting heat generation and efficiency 
as a function of the error and the erasure time. We stud- 
ied the behavior of the erasing power which characterizes 
the average erased information per unit time achieved by 
generating a given amount of heat. We finally studied 
the erasing efficiency at maximum erasing power. 

In the second part of the paper we considered infor- 
mation erasure with help of a feedback process. In this 
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